Room-temperature implementation of the Deutsch-Jozsa algorithm with a single 

electronic spin in diamond 
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The nitrogen- vacancy defect center (NV center) is a promising candidate for quantum information 
processing due to the possibility of coherent manipulation of individual spins in the absence of 
the cryogenic requirement. We report a room-temperature implementation of the Deutsch-Jozsa 
algorithm by encoding both a qubit and an auxiliary state in the electron spin of a single NV center. 
By thus exploiting the specific S = 1 character of the spin system, we demonstrate how even scarce 
quantum resources can be used for test-bed experiments on the way towards a large-scale quantum 
computing architecture. 
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Quantum computing (QC) outperforms its classical 
counterpart by exploiting quantum phenomena, such as 
superposition of states, entanglement and so on. Al- 
though the rudiments of QC are clear and some quantum 
algorithms have been proposed so far, implementation of 
QC is still experimentally challenging due to the decoher- 
ence induced by coupling to the environment. To avoid or 
suppress the decoherence, operations on most QC candi- 
date systems that are considered scalable to a large num- 
ber of qubits have been carried out at low temperatures. 
Despite that technical effort, only a few quantum gate 
operations could be achieved coherently within a single 
implementation. 

Compared to other QC candidate systems, however, 
the nitrogen-vacancy defect center (NV center) in dia- 
mond is an exception where QC operations on individual 
spins could be achieved at room temperature [lj. Since 
the first report for optically detected magnetic resonance 
on single NV centers in 1997 [2|, much progress has been 
achieved in exactly manipulating this system. As both 
electronic and nuclear spins are now well controllable 
0-0], NV centers could be used as very good building 
blocks for a large-scale QC architecture. In the QC im- 
plementation, the electronic spins are manipulated in an 
optical fashion, and the nuclear spins are operated by 
means of hyperfine coupling. Currently available tech- 
niques have achieved the quantum information storage 
and retrieval between electron spin and the nuclear spins 
0| . This technique also enables rapid, high-fidelity read- 
out of quantum information from the electron spin 0]. 
On the other hand, using the nuclear spins and additional 
electron spins as a controllable environment, a 'surpris- 
ingly different behavior' in the dynamics of the single 
electron spin was observed in different situations [l(| • As 
a result, NV centers are considered as an excellent test 
bed for models and schemes of QC. 



Despite all this progress in quantum gate realization, 
no real quantum algorithms have yet been demonstrated. 
In this Letter, we report a room-temperature implemen- 
tation of a quantum algorithm, i.e., the refined Deutsch- 
Jozsa (RDJ) algorithm (l7j . using a single NV center. 
The RDJ algorithm is the simplified version of the orig- 
inal DJ algorithm [ll|, one of the most frequently men- 
tioned quantum algorithms. As the first proposed quan- 
tum algorithm, the DJ algorithm has been employed in 
different systems to demonstrate the exponential speed- 
up in distinguishing constant from balanced functions 
with respect to the corresponding classical algorithm. 
For example, it has been carried out experimentally in 
nuclear mag netic resonance systems fl2| , in quantum dot 
systems liJLLJ], by linear optics 15j, and by trapped ions 
16]. Compared to the original DJ algorithm, the refined 
version [171 ] removes the qubit for the evaluation of the 
function, which remains unchanged during the algorithm 
implementation. As a result, it can reduce the required 
qubit resources, but still maintain the superiority due to 
quantum power over the corresponding classical means. 

We realize the RDJ algorithm by encoding both qubit 
and an auxiliary state in the S = 1 electron spin of an 
NV center. To the best of our knowledge, this is the 
first room-temperature implementation of a quantum al- 
gorithm on individual spins. To carry out the single- 
qubit RDJ, we need Hadamard gates and f-controlled 
gates. The former produces the superposition of states 
from the input state and, after the evaluation function 
has run, reconverts the superposition to a detectable po- 
larization output. The f-controlled gate is defined as 
Vf\z) = (— l)-^ z )|z), where z =0, 1, and f(z) is embod- 
ied by four functions with fi(z) = and f%{z) = 1 for 
constant functions, and with fa(z) = z and fi{z) = 1 — z 
corresponding to balanced functions. As a result, for a 
two-level system, Vf i can be written explicitly as Vf t — 
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FIG. 1: (a) Atomic structure of the NV center in diamond [J. 
The NV center comprises a substitutional nitrogen (N), and 
a neighboring vacancy (V). (b) Energy level diagram of the 
electronic ground state showing the zero field splitting, (c) 
Fluorescence microscopy image of the single NV center, (d) 
Optically detected magnetic resonance spectrum for the single 
NV center obtained by a frequency sweep using 4000 averages. 
Lorentzian peaks (red line) were fitted to the experimental 
spectrum (black circles). 

-V h = (J?) and V f3 = -V fi = (5_°i). For the two 
levels |0) and |1), the qubit is initially prepared in |0). 
After a Hadamard gate, followed by the f-controlled gate 
Vf i , the state of the system evolves to 

-^[(-i)^(°)|o) + (-i)^( 1 )|i)] ! 

with j =1, 2, 3, 4. So the constant (balanced) func- 
tion, after the second Hadamard gate, will evolve to 
( — |0(1)), which can be identified by a single read- 
out. To simplify the experimental realization, we replace 
the Hadamard gate by a selective 5 microwave pulse, 
which inverts the output of the refined DJ version com- 
pared to the original algorithm scheme but otherwise 
does not change our conclusions. 

The structure and the ground state of the NV center we 
employ are depicted in Fig. [TJ where the defect includes a 
substitutional nitrogen atom and a vacancy in the nearest 
neighbor lattice position (Fig. Ufa)). It is negatively 
charged since the center comprises six electrons, two of 
which are unpaired. Our sample is a commercial diamond 
nanocrystal (nominal diameter 25 nm). Fig. QJc) shows 
an image of the nanocrystal detected by the fluorescence 
microscopy. Since we did not apply an external magnetic 
field, the Hamiltonian of the NV center is given by [20| : 

H = sBs + S*£t, (1) 

where S and I are the spin operators associated with 
electron and nucleus, respectively. The second term rep- 
resents the hyperfine interaction between these spins, 
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FIG. 2: Transient nutation of the electron spin between 
ground state sublevels of the single NV center, where the up- 
per plot is for nutation between m s — and m 3 — — 1 sub- 
levels and the lower one for the nutation between the m a — 
and m 3 — 1 sublevels. The experimental data (black solid 
circles) were fitted by a damped sine function (red line) writ- 
ten as y — yo + Ae~ Rt cos(cot), with the intensity offset yo and 
the amplitude A to be 0.5. This normalization is independent 
from the decay rate R and nutation frequency cj. The fast 
Fourier transform (FFT) was applied to obtain the Rabi fre- 
quencies, (a) Nutation experiment at frequency 2.8254 GHz 
(MWi) with 1 million averages. Its FFT spectrum (b) shows 
a 7.87 MHz Rabi frequency, (c) Nutation experiment at fre- 
quency 2.8644 GHz (MW 2 ) with 1 million averages. Its FFT 
spectrum (d) shows a 4.26 MHz Rabi frequency. 

which is not employed in our operations but is the 
source of electron-spin dephasing. The optically de- 
tected magnetic resonance spectrum of the NV center 
(Fig. Hid)) shows that the magnetic hyperfine interac- 
tion (A ~ 2 MHz [IH) was not resolved in our experi- 
ments. The first term is the zero-field splitting (or fine 
structure) caused by mutual interaction of the two un- 
coupled electrons. This dipolar term can be written as 
H D = D[S 2 Z - lS(S + 1)] + E{S 2 X - Si). For diamond 
nanocrystals, the value of E is usually non-zero because 
of the strain induced by its vicinity to the surface, which 
assures that all degeneracies of the triplet ground state 
are lifted. We find D = 2.8449 GHz and E = 19.5 MHz 
(Fig. lib)), corresponding to a splitting of 2.8644 GHz 
between the states |0) (m s = 0) and |1) and 2.8254 GHz 
between |0) and | — 1). In our experiment, we encoded 
the qubit in |0) and | — 1), and took the level |1) as an 
auxiliary state. 

The experiments were carried out with a home-built 
confocal microscope operated at room temperature. The 
sample, mounted at the focus of the microscope, is il- 
luminated by a diode-pumped solid-state laser (Oxxius, 
SLIM-532S-50-COL-PP) at a wavelength of A = 532 nm. 
A piezoelectric scanner (Physik Instrumente, P-562.3CD) 
was used to control the focus of an oil immersion objec- 
tive (Olympus, PlanApoN, 60x, NA= 1.42). The NV 
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center fluorescence is separated from the excitation laser 
with a long wave pass filter (Semrock, BLP0f-635R) and 
then collected by a silicon avalanche photodiode (APD) 
(Perkin Elmer, SPCM-AQRH-f 3). We constructed two 
synchronized microwave channels that provide the setup 
with the ability to output microwave pulses with two dif- 
ferent frequencies. The microwave is coupled to the sam- 
ple by a 20 /im diameter copper wire acting as antenna. 
The whole system is orchestrated by a word generator 
(SpinCore Technologies, PBESR-PRO-350). 

As a preparation for the RDJ algorithm, we have 
first accomplished coherent spin resonance between the 
ground state sublevels. Fig. [2] shows the transient nu- 
tations of a single NV center. The initialization of state 
|0) with > 90% probability is achieved by a 5 us green 
excitation, followed by a waiting time of 5 fxs [19J. A 
microwave pulse of variable duration was then applied to 
the NV center, and the spin state was read out by moni- 
toring the fluorescence intensity. The experimental data 
show a periodic modulation of the fluorescence signals 
between the |0) and | — 1) (Fig. Ufa)) and between the 
|0) and |1) (Fig. Hfb)), respectively. From the figures, 
we can extract the Rabi frequencies under microwave ir- 
radiation ((Fig. [5Jc) and (d)), which is important for 
performing the gates in the RDJ algorithm below. We 
used 64 ns and 118 ns long tt pulses for the two mi- 
crowave channels, respectively. The decay of coherent 
oscillations is due to electron spin dephasing, affecting 
the visibility of the spin state in a predictable manner. 

In the implementation of RDJ, after preparing the ini- 
tial state |0), we applied a selective J microwave pulse 
in the MW\ channel, yielding the superposition. The 
following f-controlled gate operations have been imple- 
mented by combinations of 2ir pulses in the four possible 
cases (Fig. [3]). Making use of the auxiliary state |1), we 
applied a 2ir pulse in the MW2 channel, introducing a 
tt phase shift on the state |0). This is equivalent to a 
7r rotation about Z axis in the subspace spanned by |0) 
and I — 1) [21] . Switching the 532 nm laser on again, the 
final result was read out from the fluorescence collected 
by the APD. Fig. [3] shows the four groups of microwave 
pulse combinations we used, corresponding to the four f- 
controlled gate operations. The intervals between pulses 
are set to zero in our experiments and we merged some 
pulses into a longer one in order to simplify the pulse 
sequence. For example, a MW\ tt pulse was used in 
Figj3^a) instead of the two ? pulses. This increases the 
selectivity of these pulses, but the pulses are still broad 
enough to cover the resonance line. 

The results of the RDJ algorithm are shown in FigJU 
Each point (gray stars) represents an individual experi- 
ment with 50 million averages, and we did four sets of 
experiments corresponding to the four f-controlled gate 
operations. Note that since we have used selective § 
pulses in place of the Hadamard gates, the output is in- 
verted with respect to the original algorithm. The weak- 
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FIG. 3: Diagram of the experimental pulse sequences used to 
realize the RDJ algorithm. The 532 nm laser (green line) was 
used to initialize the state of the NV center to |0) and was 
shut off during the algorithm. Later, the laser was switched 
on again for detection. MWi (red line) and MW2 (blue line) 
are two microwave channels which excite different transitions 
selectively. The first MW\ ^ pulse was used to generate a 
superposition state in the qubit and the last MWi \ pulse was 
taken for detection by reconverting coherence to population 
differences. The V/, operations (i = 1, 2, 3, 4) were realized by 
combinations of (both MW\ and MW2) 2n pulses, as shown 
between the vertical dashed lines. 
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FIG. 4: The output of the RDJ algorithm is illustrated by 
the intensity of the fluorescence, where (a) and (b) with the 
weakest intensity correspond to the constant function, and 
(c) and (d) with the strongest fluorescence indicate the bal- 
anced functions. While the original data (gray stars) suffer 
from dephasing during the pulse implementation, the outputs 
are clear enough to demonstrate the differences of f-control 
functions. The black circles are plotted with the dephasing 
effect in spin nutation compensated. 
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est fluorescence intensity (minimum population in |0)) 
thus corresponds to the constant functions V/i and Vf2, 
and the strongest intensity (maximum population in |0)) 
indicates balanced functions V/3 and Vf^, The fluores- 
cence difference obtained for constant and balanced func- 
tions is 56.9% (gray stars in Fig. 2]), which is reduced 
mostly due to dephasing during the pulse operations (es- 
timated as ~ 59.6% from Fig. [2]). Nevertheless, it is clear 
enough to illustrate the accomplishment of the RDJ. For 
clarity, we have compensated the dephasing effect using 
the results in Fig. [5] As a result, a small deviation of the 
compensated experimental data (black circles in Fig. [4j 
from the theoretical expectation (red lines in Fig. |4} re- 
mains, due to operational imperfections of the microwave 
pulses. 

In contrast to room-temperature experiments with nu- 
clear magnetic resonance using spin ensembles 
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our 

experiment works on a solid-state quantum system. As 
a result, we have achieved a pure-state QC implementa- 
tion at room temperature. With respect to other sys- 
tems 0, m for coherently manipulating individ- 
ual spins, our implementation without cryogenic require- 
ments greatly reduces the experimental challenge for car- 
rying out QC. Moreover, as our qubit in the NV center 
can be fixed and manipulated exactly, the QC operation 
in our case is deterministic and efficient. 

The full demonstration of the power of quantum algo- 
rithms requires large-scale QC. Optical coupling of spa- 
tially separate NV centers might be achieved by putting 
the centers in optical cavities, which enhances both the 
zero phonon line and the collection efficiency of the emit- 
ted photons. Considerable efforts have been made for 
fabricating thin, single- crystal diamond membranes [22| , 
whispering-gallery mode resonators [23| and photonic 
band-gap microcavities 24]. However, since these sys- 
tems usually work well only at low temperatures, and 
since the excited states of the NV centers are not well 
protected from decoherence at room temperature due to 
spin-orbit coupling [25j , cryogenic operation seems neces- 
sary for extending the NV-center QC architecture in this 
way. Therefore, how to accomplish a large-scale QC at 
room temperature is still an open question [l[ . Moreover, 
once more qubits are involved in the system, the required 
operations get more complicated and time-consuming. 
This implies that we need to fasten the operations or to 
effectively suppress decoherence. Some first explorations 
into these aspects have been reported 0, E3| ■ Neverthe- 
less, our present experiment has clearly shown the unique 
opportunity provided by NV centers to study the physics 
and application of single spins, and also demonstrated the 
great potential of the NV system for QC. 

In summary, we have accomplished a RDJ quantum al- 
gorithm using only the electron-spin of a single NV center 
at room temeprature by exploiting the S = 1 character 
of this system. Our experimental data (after compensa- 
tion of a systematic dephasing effect) fit the theoretical 



prediction well with a small deviation due to pulse im- 
perfections. Although building a scalable quantum com- 
puter is pretty hard with current technology, successful 
implementations of existing quantum algorithms using 
available QC building blocks would be definitely helpful 
in stimulating inventive new ideas and further technolo- 
gies. In this sense, it is of great importance for our ex- 
periment to demonstrate the power of QC at room tem- 
perature in a real solid system, using minimal quantum 
resources. The possibility of carrying out quantum su- 
perposition and interference at room temperature makes 
future work toward large-scale room-temperature QC ar- 
chitectures worthwhile. 
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